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1.0 


INTRODUCTION 


The purpose of the Shuttle LOX Loading Transient Study is to provide 
a technical basis for configuration and operation of the LOX loading 
facilities to be used for the Space Shuttle Main Propulsion Test Articl 
(MPTA) at the National Space and Technology Laboratories (NSTL), and 
for the Space Shuttle launch facilities at the John F. Kennedy Space 
Center (KSC). The principal objectives of this study are to determine 
the acceptability of these LOX loading facilities for each proposed 
application and to determine the degree of similarity between the 
systems with respect to the applicability of data developed during 
operations at NSTL to the launch site operations at KSC. Further 
objectives are to determine facility modifications which would improve 
the similarity. 

This study was divided into three major tasks: 

(1) Cryogen Transfer Computer Program Development and 
Verification, Reference 1. 

(2) National Space and Technology Laboratories (NSTL) LOX 
Loading Facility Analyses, Reference 2. 

(3) Kennedy Space Center (KSC) LOX Loading Facility Analyses, 
Reference 3. 

During the course of the NSTL and KSC facility analyses, alternate 
facility configurations and operations were considered and evaluated. 
This report presents a summary of the overall Shuttle LOX Loading 
Transient Study results, final data comparisons, conclusions, and 
recommendations. 

Contract NAS8-30745 was modified for an additional major task to 
investigate the feasibility of Transient Cryogen Transfer Program" 
(TCTP) application to the Space Shuttle Main Engine (SSME) hydrogen 
system analysis. The results and conclusions from this task are 
being published in a separate report. 
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2.0 COMPUTER PROGRAM DEVELOPMENT AND VERIFICATION 

The Shuttle LOX Loading Transient Study was based upon further develop- 
ment and verification of the "Transient Cryogen Transfer Program" (TCTP) 
and "Steady State Cryogen Transfer Program" (SSCTP) as documented in 
Reference 1. During the NSTL and KSC facility analyses, extensive TCTP 
modification and development was required to simulate the proposed 
system configurations and operations. These modifications were incorporated 
into the up-dated program user's manual presented in Appendix A. The 
TCiP modifications provide a more generalized capability to perform 
transient analyses of various cryogen transfer system configurations 
and operations with various known boundary conditions. 

During the NSTL facility analyses, the TCTP was improved to include the 
following capabilities; 

(1) Throttlable flow control valves. 

(2) Storage tank bottom pressure as a function of storage tank 
ullage pressure, liquid level, and suction line flow rate. 

(3) Pump inlet pressure and NPSH as a function o\ storage tank 
bottom pressure and suction line flowrate. 

(4) Pump by-pass flow as a function of pump discharge and storage 
tank pressures. 

(5) Tank head flow simulation starting from zero flow rate into 
hot (ambient temperature) lines. 

These TCTP improvements are discussed in Reference 2. The tank head 
flow capability was verified by analyzing the Skylab S-IB Stage LOX 
loading system initial tank head flow chi 11 -down and comparing 
analytical results to available measured data. A comparison of 
analytical and existing measured data is shown in the Appendix B 
presentation. Considering the accuracy that the available strip 
chart data can be read (jL 2 psi and +^1.5 LB/SEC) and the spread 
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2.0 (Continued) 

between the two redundant pump discharge flow rate transducers (33G01 
and 33G08) the analytical results agree well with measured data. This 
analysis also revealed that transient quality and flow rate gradients 
exist in the system at the end of the tank head flow portion of line 
chi 11 -down just prior to the initial 1 ,000 GPM pump RPM ramp. Using 
the transient flow rate, quality, pressure, and temperature gradients 
at the end of tank head flow as input, the pump flow portion of the system 
chill-down was reanalyzed. Results of this additional program verifi- 
cation are shown in Appendix B. The resulting transient flow rate and 
pressure gradients were approximately the same as those reported in 
Reference 1. The analytical time at which the vapor volume fraction 
goes to zero at the LOX mast measurement location is 212 seconds from 
the initial RPM ramp, while the measured data indicates that the transition 
to liquid flow occurs at approximately 203 to 224 seconds, as reported 
in Reference 1. 

During the KSC facility analyses. Reference 3, additional improvements 
were made to the TCTP to include the following capabilities: 

(1) Throttlable branch flow to the supply tank reservoir (KSC 
LOX by-pass). 

(2) Check valve downstream of pump chill-down drain line. 

(3) Improved system bleed flow simulation to compute bleed flow 
rate as a function of system node pressure and density and 
bleed system exit pressure. 

(4) Option to compute required pump speed as a function of input 
system entrance flow rate and computed pump discharge 
pressure instead of system entrance flow rate as a function 
of input pump speed and computed pump discharge pressure. 


2-2 


D180-19190-1 
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3.0 SUMMARY AND COMPARISON OF RESULTS 

The suggested tank head flow facility chi 11 -down period is 4 minutes 
for the NSTL facility, Reference 2, and 10 minutes for the KSC 
facility, Reference 3. The chi 11 -down times are expected to be 
different for the two facilities due to differences in system configu- 
rations, system operations, heat loads, and supply tank conditions. 
However, the transient conditions in the Orbiter and 17 inch ET fill 
duct are very simular over the facility chi 11 -down period. From 
Reference 2, Figure 2-3, and Reference 3, Figure 3-2, the Orbiter 
(ORB)ZGround Service Equipment (GSE) interface flow rates are between 
2 and 4 LBn^/SEC over most of the chili -down period at both facilities. 
At the end of the tank head flow chi 11 -down period, the ORB/GSE 

flow rates are 3.6 and 2.9 LBn^/SEC for the NSTL and KSC facilities. 

During this period the maximum, or peak, gas flow rate to the ET 

is 3.0 and 8.5 LB|^/SEC at NSTL and KSC respectively. Comparison of 

Reference 2, Figure 2-4 and Reference 3, Figure 3-3 shows that the per- 
cent vapor by mass (quality) decreases slightly more at NSTL than at 
KSC during the facility chi 11 -down. At the end of this period the 
NSTL ORB/GSE quality is 0.14 and the KSC ORB/GSE quality is 0.26. From 
Reference 2, Figures 2-5 and 2-6, and Reference 3, Figures 3-4 and 3-5, 
the ORB/GSE interface and vehicle fluid and duct wall temperatures are 
essentially the same for both facilities at the end of tank head flow 
chi 11 -down. 

Reference 2, Figure 2-9, shows that the NSTL facility is filled with 
liquid with 100 percent liquid flow at the ORB/GSE interface at 1.9 
minutes after the pump flow slow fill initiation and the 17 inch 
duct is filled to the upper aluminum/stainless steel joint at 16.5 
minutes after slow fill initiation. Reference 3, Figures 4-4 and 
4-11, show that the KSC facility is filled with liquid at 2.8 minutes 
after slow fill initiation and the 17 inch duct is filled to the 
upper aluminum/stainless steel joint after a total of 24.7 minutes 
of slow fill time. The longer required fill time for the KSC facility 
is due to the higher overall heat load for the facility. Alternate 
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3.0 (Continued) 

KSC Orbiter chill down and vehicle slow fill analyses with various 
flow rates showed that the required fill times may be changed 
significantly with different controlled flow rates. At the end of 
the slow fill periods evaluated, Reference 2, Figures 2-10 through 2-12, 
and Reference 3, Figures 4-12 through 4-14 show that the ORB/GSE 
interface and vehicle pressures, fluid temperatures, and duct wall 
temperatures are all very similar. 

Results of these analyses show that the facility and Orbiter chill - 
down and slow fill can be accomplished within the target time of 
43 minutes at both facilities provided that the SSME engine chill- 
down can be accomplished simultaneously with the Orbiter chi 11 -down 
and vehicle slow fill, and no engine chill hold time is required 
at KSC, 
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4.0 CONCLUSIONS AND RECOMMENDATIONS 

The results of the NSTL and KSC facility analyses, References 2 and 3, 
lead to the following major conclusions: 

(1) Shuttle vehicle interface conditions at both facilities will 
be approximately the same after the tank head flow facility 
chi 11 -down period. 

(2) Similar vehicle interface conditions can be maintained at 
both facilities during and after the pump flow slow fill 
initiation with flow control valve or pump RPM adjustments 
at either facility. 

(3) The facility/vehicle chill-down and slow fill to the ET 
(tank bottom) can be accomplished within the target time 
of 43 minutes at both facilities. 

(4) Desired flow control or operational procedure adjustments 
and accurate system performance predictions for KSC launch 
operations can be derived from NSTL test results coupled 
with further TCTP analyses. 

Results of the NSTL facility analyses, Reference 2, indicate that a 
sudden increase in flow control valve upstream and pump discharge 
pressures will occur when 100 percent liquid flow reaches the flow 
control valve (LCV-114). At this time, the flow downstream of LCV-114 
will be choked. This sudden increase in pressures and choked flow 
condition could be avoided by one or more of the following changes to 
the planned system configuration or operation: 

(1) Increase the barge tank ullage pressure, system flowrate, and 
tank head flow period to achieve 100 percent liquid at LCV- 
114 prior to pump start. 

(2) Increase the replenish pump RPM slowly with LCV-114 held 
open until 100 percent liquid reaches the valve at a low 
flowrate. 


4-1 


D180-19190-1 


4.0 (Continued) 

(3) Use a larger flow control valve and downstream line size. 

Results of both facility analyses indicate that after the 17 inch 
Shuttle ET fill duct is filled with liquid to the aluminum/stainless 
steel joint (10 feet below the ET) steady state conditions will be 
approached with the liquid mass flow rate approaching the total mass 
flow rate (low quality) at the ET entrance. However, the 2-phase 
fluid will be approximately 60 percent vapor by volume at the ET 
entrance. If it is desired to guarantee 100 percent liquid flow 
at the ET entrance at the end of the system chill-down and slow fill, 
one or more of the following changes should be considered for the 
planned system configuration or operation: 

(1) Decrease the system inlet or vehicle interface LOX 

temperature. 

(2) Decrease the external heat load (improve the system insulation). 

(3) Increase the system flow rate during the 17 inch duct fill. 

(4) Maintain a higher ullage pressure (approximately 19.5 PSIA) 
in the ET to increase the tank bottom pressure above 
saturation pressure. 
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APPENDIX A 

ENGINEERING DESCRIPTION AND USERS GUIDE 
FOR THE CRYOGEN TRANSFER PROGRAfYS 

This Appendix presents an up-date to Boeing Letter Report 5-9030-HT-1 
"Cryogen Transfer Computer Program Development and Verification", 
Appendix B, dated September 3, 1974. Contents of this Appendix 
represent the latest Cryogen Transfer Program theory and input data 
format for both the TCTP and SSCTP. 
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ATTACHMENT 1 


TRANSIENT CRYOGEN TRANSFER SYSTEM ANALYSIS 

1.0 CONTROLLING EOUATIONS 

1.1 f1ASS 

From Rfiference 1, the equation of continuit'' for nt.‘*-riirif'n';icmHl trnn 
sient flow is, 


A- 

3t 




or in finite difference notation for constant flow rates over ,‘t, 



At 


! 



( 1 ) 


where at any instant the mass contained within the control volume or 
line node is, 




( 2 ) 


and the flov/ rate leaving a node during a time step is, 

Mf.| = 

where p-'m is the effective average density over the time step, and for 
the finite difference integration, assumed to equal the final value for 
the time step. 

Pn PFN 


(3) 


( 4 ) 


1.2 ENERGY 

The thermodynamic system boundary prior to a time step, dt, is considered 
to include the fluid within the control volume, or line node, plus the 
increment of fluid v/hich will enter the node during the time^step. The 
same system at the end of the time step includes the fluid within the 
line node plus the increment of fluid which left the node during the 
time step. In the limit as dt approaches zero, the system boundary 
coincides with the control volume, or line node. 


The energy added to the system during the process is, 

where (P,^ -,3) is the flow work done on the fluid within the node 

per pound of fluid forced past the entrance, and puj is the rate of heat 
transfer from the v/all to the fluid. 
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1.2 (Continued) 

The energy removed from the system during the process is the flov/ work 
done by the fluid within the node to force Mj>|dt pounds of fluid past 
the node exit. 


The energy possessed by the system in its initial state is, 


g J 2g 

L» 


and in its final state, 

MpUp + (f\^dt)(U + 

Equating the difference between the energy added and removed to the change 
in energy possessed by the system gives, 

? 2 

= 1^1 U _ M II + Cm CII 4- SA- 4- 1 \ - (M Cl! 4- 3A^ + ^ — 


pup 


U. + (M.,dt) (U + 


gcJ 2g^J 


)n - ^ M 


By definition the enthalpy is, 
h MU + 

and the final mass from equation (1) is: 

Mp = (f^N-l ■ %' 

By combining terms and substituting the above expressions for enthalpy and 
final mass, the energy equation may be solved in finite difference form 
for the final internal energy, 2 2 

U lg^^N-1 " ^ ^g^'^N] At 

[Mj + - l\) At] 

^ (hcN ^NW (Tw - T)n At/36Q0)+(MU),i 
[Mi + (l^|,,_i - M|^P At] 
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1.2 

v/here 

1.3 

1.3.1 


(Continued) 

" 3^00 

MOMENTUM 
Constant Area 


For a constant area control volume or line segment and one-dimensional 
transient flow, the momentum theorem, Reference 1, becomes 

>'Fm " - <^N-1 

where the time rate of change in momentum for node, n, is evaluated at 
the node entrance conditions. The partial derivative of (pV) with respect 
to time is. 


av 


it (pv) = p S + 

where p and V are instantaneous local values. For finite time steps, the 
partial derivative of pV with respect to time becomes. 


Pc^Pt 1 

[ 2 (vp - Vj) + ■ 2 ^^F " At 


= (ppVp - Pt^t) 


1 


VI' At 


The line element or node volume is defined as, 

V|y| = 

The momentum theorem may then be written in finite difference form as, 

^^N " i^^^CN-DF^VllF “ P(N-1)l''(N-1)I^^N*^N 

2 „ 

+ (*^N'^N ) A|^ At - (p|\|_i 

The summation of forces in the direction of flow on a constant area line 
increment is. 


aPA 


V/h>F - 'n-/n ^ -3X "’Pp ^ ■''“n 

* *‘’NF*gr' 


dl 


N 
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1.3.1 (Continued) 

For a finite line node length, the shear force ( ! is balanced 

by the pressure drop due to friction along the constant area nine wall, 
i . e . , 


(TttD^jAS^P - = 0 

O'" ^PpN ' ^'"n-I ■ 

where the subscript, F, denotes frictional nressure drop. 


The summation of forces over a finite line length may now be written, 

>-Fn “ -''n [<Pn - PS-l) " ''PfN " - Vl>l 


By combining terms and rearranging, the momentum of equation may now be 
solved for the change in pressure across a constant area line length. 


(Pn " ^N-1^ " ^P(N-l)l''(N-l)I ■ P(N-1)F'^(N-1)F^ 


- ^^PfN - (Z, - Z,_^) (6) 

1,3.2 Area Change 

If a constant area node is preceded by an area change, the velocity and 
pressure entering the node must be corrected. It is assumed that the mass 
flow rate entering the area change interface is equal to the mass flow 
leaving, and that negligible change in density occurs across the thin 
interface. The flow conditions across an area change may be approximated 
by the Bernouli equation modified to account for frictional pressure 
drop due to abrupt area changes, 

^''n-1 “ '"n-1^ ^ Pn-1 ( ^N-1 - ^N-1 ) = 


144 


The velocity entering node, n, is then 


A- 5 


''N-1 "''N-1 


(7) 
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1.3.2 


(Continued) 


and by substitution into the modified Bennouli equation, the clianqe in 
pressure across an area change just upstream of node, n, is 




- 


1-1 


) = [1 - 


C ‘ - 1 




( 8 ) 


144 
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TRANSIENT CRYOGEN TRANSFER SYSTEM ANALYSIS 

2.0 SUBSIDIARY EQUATIONS 

2.1 HEAT TRANSFER 

2.1.1 Heat Balance 


The heat transfer rate from the external environment to the wall minus 
the heat transfer rate from the wall to the fluid is equal to the rate 
heat is stored by the wall: 


^EXT 


- q = 


At 




Solving for the final node wall temperature, 

T - T + ^CN^NW^’’'WNI ■ ’’’nI^ 

TwNF - %I ^ M^C^ - 3600 (9) 

The transient heat transfer rate q^vj is determined as a function of 
the initial wall temperature for ^ the time step by the following 
equation: 

*^EXTj^ " ^EXT(MAX)f^ (520-T^^j )/(520-T^ p) (9a) 

where qryr/MAY^ ■'S fhe maximum steady state external heat transfer 
rate Tiquid in the line and with the wall 

temperature, T^, approximately equal to the pump discharge temperature, 

^ is input to the computer program. 


2.1.2 Two Phase Flow Heat Transfer Coefficient 

For this problem, heat transfer from the wall may take place to a single 
phase (liquid or gas) or a two phase fluid. Also, the quality of the 
fluid in a line mode may vary with time. Since the available tv^o phase 
correlations do not reduce to the single phase liquid correlation as 
the quality approaches zero, a linearization between the two correlations 
will be used. The selected two phase correlation does, however, roughly 
approximate the single phase correlation at high quality. 

The two phase heat transfer correlation from Reference 2 was selected as 
the correlation most applicable to the cryogen transfer system with 
initial chi 11 -down and significant variations in hydrostatic head which 
cause liquid vaporization. Although the system being analyzed will 
have relatively low wall to fluid temperature differences, it is reasoned 
that the two phase flow model will be similar to that envisioned by the 
authors of Reference 2 due to the relatively high liquid velocity at 
significant vapor volume fractions. The flow model is assumed to have 
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2.1.2 (Continued) 


a vapor film at the wall with liquid droplet impingement on the wall 
at higher values of vapor volume fraction. The following Nusselt number 

two-phase flow at qualities between 0.010 and 1.0 
was obtained from Reference 2 and a telephone conversation with 
nr. K. c. Hendricks correcting a printing error in Reference 2: 


NU^ 

^^"Oll + 1.93'i'^.j. 


1 


+ 0.15* 


( 10 ) 


where is 


¥ 


TT 


the Martinelli two-phase flow parameter, 
- /'Vv 0.5 0.1 


( 11 ) 


*The above constant (0.15) was changed to 0.221 to improve correlation 
with S-IB Stage LOX loading system two-phase heat transfer. 


The reference or calculated Nusselt number used in the Reference 2 
correlatijin is evaluated with a Dittus-Boelter equation using the two 
phase bulk density and the vapor film viscosity, specific heat, and con- 
ductivity For qualities between 0.010 and 1.0, the two-phase flSw 
film coefficient is then given by: 


'CN 


'CLVN 


0.021 


Ku (pVD) 


D 


[- 


>‘V 


Nj 0.8 


(Pd ) 


0.4 


(!?!£) 

'uc ' 


( 12 ) 


using average vapor film properties 
Additional discussions of this two- 
are given in References 3 and 4. 


over the expected temperature range, 
phase flow heat transfer correlation 


For = 0, the single phase liquid film coefficient is given by: 

0 

RL‘ 


h = h = n noi 0 /^L N N\0.8/n \0.4 

TN \ln ^ ^ ) 


1 


( 13 ) 


using the average liquid properties over the expected temperature range, 


For X[^ between 0 and 0.010, a straight line interpolation between the 
single phase liquid film coefficient and the 2-phase value is assumed 
to be approximately correct. 
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2.2 FRICTIONAL PRESSURE DROP 

2.2.1 Two-Phase Flow Line Losses 


The selected two-phase flow line frictional pressure drop correlation was 
taken from Reference 5 which uses a two-phase flow friction factor 
correlation from Reforonce 6. No-slip homogeneous bubble or froth flow 
is assumed due to the relatively high anticipated liquid velocities. 

For no-slip flow, the vapor volume fraction is: 


1 + ^ (^) 

Pl 


The two=phase fluid viscosity is: 


(14) 


^N " '"V “N ^L " “N^ 


where average values of liquid and vapor viscosities ma'/ be used over the 
expected temperature range. As suggested by Reference 5, a single 
phase friction factor is determined from an existing single phase 
correlation as a function of the two-phase Reynolds' number : 


fQ,^ = 0.0014 + - 


PVg 0.32 
^ u ^N 


(15) 


Tin's single phase correlation is for smooth pipes. Other existing 
correlations may be substituted for other than smooth pipes. The tv/o- 
phase flow friction factor correction from Reference 6 is: 


/tPx ^ /1/2 EuTP ^ 

^f ^N ^ \ f ) - 

0 ^0 

P (l-«) * 

^ “ 1 .281 + .4781n(l-a) + .4441n2(l-a)+.0941n3(l-a) + .008431n4(l-a)^ 

N 


The friction pressure loss in each constant area node is: 


AP 


FN 


= (hi) TK + 4 foil 


/^N-lVN?h 
^ 288 g. ^N 


(17) 


where is the loss coefficients for hardware elements in each node. 

*Multiplying (fjp/f^)-l by a factor of 3.25 resulted in a slight improve- 
ment of analytical correlation with measured S-IB Stage LOX loading system 
two-phase flow pressures. 
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2.2,2 Line Area Change Losses 

Reference 7 suggests that for homogeneous two-phase flow, head losses 
across abrupt expansions and contractions can be predicted using co- 
efficients for single phase flow. For contractions: 


AP 


, _ *^AN /^N-1''N-1n 

FN T?? ^ 2g^^ 


(IH) 


From Reference 8, the contraction loss coefficient as a function of 
diameter ratio is; 

TABLE 2-1 


A 

'n-1 


""AN 


0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

o.a 

0.9 

0.45 

0.42 

0.39 

0,36 

0.33 

0.28 

0.22 

0.15 

0.06 


For Expansions 

"Pfn 


1^ 2 
^ AN ^fN-l^N-l 


T44 




4 


(19) 


where the expansion loss coefficient from Reference 8 is; 

K - h 
^AN “ 


2.3 


"N 

SYSTEM PRESSURES 


( 20 ) 


The total change in pressure across each node and an area change just 
upstream may be determined by adding equations (5) and (7): 


AP 


N 


^*^N ■ '^N-1^ ^ ■ '^N-1^ ^’’n-1 “ ^N-1^ 


(24) 


The pressure leaving each node, n, is then 

K 

P,^ = P|< + j: (-AP), 


N+1 




where P|^ is the known final node discharge pressure. 


( 25 ) 
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2.4 SYSTEM FLOW RATES 

2.4.1 System Entrance Flow Rates 


In addition to the above flow controlling equations, a relationship for 
flow rate entering the system is required to complete the solution for 
system flow rates. This analysis considers the incoming flow to be 
controlled by the flow rate characteristics for a given pump as a 
function of outlet pressure and RPM. For a pressure transfer system, a 
suitable flow rate equation for a known flow controlling valve,' orifice, 
or system entrance loss may be substituted. The pump performance character- 
istics may be normalized to give: 


1^1 Pi-Ptn 

(tT> = 

IN 


(26) 


The pump flow rate is then: 
= N (M^/N) 


(27) 


where N is the known pump RPM as a function of time. The flow rate leaving 
each node is determined by a simultaneous solution of the energy, mass and 
fluid property equations. 


2.4.2 System Branch or Bleed Exit Flow Rate 

For special problems where branches or bleed flow leaves a given 
system node as a known function of the system node exit pressure, 
the bleed flow rate may be subtracted from the system node exit flow 
rate. The bleed flow rate may be calculated by the following equation: 

rtl ( 288 g^)Pn''BL 

where Abl is the branch or bleed equivalent flow area, P^x is the 
branch line or bleed system exit pressure, Zex is the branch line 
or bleed system exit elevation, and Ke is the equivalent friction 
loss coefficient. 

Combining constants, this equation becomes, 
lii "^y^Pl Pn ^V*^EX“^BP2 

where Cgp.j = (288gc Ag^/Kg) 
and Cgp 2 = (^eX'^IN^^^^^ 


A-11 


D1 CO-1 91 90-1 


2.4.3 Two-Phase Choked Flov/ 

Relations were derived to check for two-phase choi t ri t'.ow and -u.j 
the flow_ rates if necessary. However, due to t.t.r- mar:y fjthf * tiotiai i. 
sensitivity problems encountered during program dovo'h/pinent irra.,ifi< 
time was available to develop an operational [jruutan i! , (i*. 

matically adjust system flow rates to prevent excoeditaj fiM;. 

at any node. The current program does limit tiode and v; it 

velocities to computed choked flow values for tho [-lirf.rso oi i-.VI;. ]. 
node pressure drops. If choked flow is exceeded at a fi r. < . 

choked flow conditions can be simulated by control ! Ina thn r;t* .m t! 
a flow control valve input as discussed in Referenno 10. 


For a unit volume, the compressibility is by dofinit! 


f !ti : 


k' 1 dV 
V iP 

The mixture volume is, 

''ll 

Therefore, for the mixture, the compressibilitv i: 


- - 7^ (— 


■) = - ^ (V,K, . V,: ,} 


The vapor volume fraction is by definition 

''v ''l 

“ " IT -^ = (1 - «) 

''m 

Substituting into the above expression for compre''<'.i[>ili?y vii ;: 

= Kl (1 - -) + Ky. 

The bulk modulus is defined as the reciprical of rtijiua i; ; ] i r 
Therefore, 


K 


1 


and. 


:l= 0 - -) 
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2.4.3 (Continued) 

From the above definitions of compressibilitv, and bulf: modulus it may 


dp 


and from Reference 1, the speed of sound is, 

9'=044) = y^gc(144) 

The choked flow rate per unit area at Mach No. equals 1 is then, 

IVJ-A* I ' 

r “ f'M" = 


thf ®xP'"ession for the mixture bulk modulus qives 

pSss'ul^e!" "" modulir 



N 


PmN 9c(H4) 

IZn + !il 

^VN 


where the node vapor isentropic bulk modulus is. 


^VN " PVN " ^V*^N 


The results of the above equations are shown in Figures 1 - 4. 

2.5 FLUID PROPERTY RELATIONS 

In order to complete the solution of the above equations, several 

additional state point property relationships or themodynamic tables 
are required. 


For the simultaneous energy and mass balance, the vapor density and 
the saturated liquid and vapor internal energy are determined from 
thermodynamic tables at the current assumed node pressure (Pp,^ = P. ). 

The node quality at the beginning or end of a time sten is determined from 
the internal energy at the current pressure. 


y _ ^N " ^LSN 
^N " U 


VSN ■ \SN 


(28) 
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2.5 (Continued) 

For the finite difference integration, the effective average quality 
for the time step is assumed equal to the final value for the time 
step, 



(29) 


The density at the beginning or end of a time step is determined from 
the quality at the current pressure, 


'"^1 '’VSN 

If the final quality is between 0 and 1 , the final temperature is the 
saturated temperature, Tcm, at the assumed final pressure. If the 
final state is subcooled (X|\|p = 0) or superheated, (Xp = 1) the final 
temperature is determined from 


where and C = Cyy for = 1 

and and C = Cp|^ for = 0 


(31) 


Also, the enthalpy leaving the node during the time step is determined 
from 


h' II X ' FN 
= ^FN 


(32) 
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TRANSIENT CRYOGEN TRANSFER SYSTEM ANALYSIS 

3.0 SYSTEM ANALYSIS 

3.1 FIRST TIME STEP INITIAL VALUES 

3.1.1 Initial Node Properties 

The initial node pressures, temperatures, and qualities must be input to 
the program. Originally, the node bulk average quality, Xn, and the 
local quality leaving the node, X'^, were input to the program. The 
current program requires only the initial node exit quality input, and 
the node bulk and exit qualities are considered to be equal. In the 
current program solution flow, the effective average quality over a 
time step is designated by X'n and is equal to Xn at initial con- 
ditions. If the node quality is between 0 and 1, the temperature 
may be input as zero and the program will determine the saturated 
node temperature from a table as a function of the input pressure. 

The initial temperature entering the first line node, T,r, is the 
known pump outlet temperature. 


The initial density, is determined for each mode from equation (30), 
where % is the known or assumed initial node bulk quality and the 
vapor density , pvsNj is evaluated from a thermodynamic table as a 
function of the node initial temperature and pressure. The liquid density, 
calculdtod from curvG fit dsts for ths spGcified cryogGn. ThG 
initial node internal energy U[^j, may be determined for each node from 
equation (<:8), where Xn is the known or assumed initial node bulk 
quality and UysN and Ulsn are the vapor and liquid internal energies of 
the specified cryogen evaluated from thermodynamic tables at the node 
saturation conditions, Pj^j and Tsnj. If the node is subcooled or super- 
neated,_the initial internal energy is determined from equation (31) 
at the input initial temperature. 

3.1.2 Initial Node Mass, Flow Rate and Velocities 

With the initial node density from equation (30) and the known system 
geometry , the initial mass in each line node may be determined from 
equation (2). The initial pump flow rate is determined from the known 
pump performance function, equation (26) and equation (27), usina the 
initial pump outlet pressure, Pu, in PSIG and the known initial' RPM. 

For initial flow at steady state, the flow rate out of each node is 
equal to the pump flow rate, 


and the initial velocity leaving each node is determined from equation (3). 

If there is an area change at a node entrance, the entering initial velocity 
IS determined from equation (7). 
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3.1.3 Initial External Heat Transfer 

If the system is initially at a non-steady condition, the initial 
external heat transfer rate is calculated with equation (9a). If 
the system is initially at a steady state condition, e.g. zero flow, 
the external heat transfer is: 

^EXT^ “ *^CNI^EXT(MAX)f^ (520-Tf^j)/[(520-T^p) h^^j + 9eXT(MAX)j^^ 

where qEXT(MAX) maximum steady state external heat transfer 

program input. ^ 


3.1.4 Initial Wall Temperature 

If the system is initially at a non-steady condition the initial node 
wall temperatures must be input to the program. The initial wall 
heat transfer coefficient is determined from equation (12), (13), or 
p4) depending on the initial quality, and at steady state the 
initial wall temperature is 

9 ext 

%I " "^NI ^ hZTKZ 


3.2 SYSTEM SOLUTION FLOW 

After the first time step initial values have been determined, the procedure 
outlined on Figure 5 is utilized to determine the pump outlet and line 
node conaitions for subsequent time steps. 

A "DO LOOP" is used for consecutive nodes, 2 through K where K is the 
system exit node, to determine the initial vapor volume fraction, 
wall film coefficient, hcN, and wall temperature Tj^pi^j. The initial 
vapor volume fraction is determined from equation (14) at the 
beginning of each time step as a function of final quality and vapor 
density from the previous time step. The wall film coefficient is assumed 
constant over the time step and based on initial flow conditions. 

If the initial node quality is zero (single phase liquid) or less than 
0.01, the heat transfer film coefficient is determined from equation 
(13). If the initial node quality is between 0.01 and 1.0, the heat 
transfer film coefficient is determined from equations (10), (11) 
and (12). 
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3,2 (Continued) 


The pump RPM (or flow rate) and exit temperature, TfIj are determined 
from input data tables as a function time. If the pump RPM is input, 
the pump flow rate is a function of pump outlet pressure, and the pump 
outlet pressure is a function of system flow rates and pressure drop, 
an iterative solution is required to determine the final pump pressure 
for each time step. For the first assumed pump pressure, the pump 
pressure rise divided by the RPM squared is assumed to be approximately 
constant from one time step to the next. Therefore, the first assumed 
pump pressure is: 


FAl 


= (P 


II 


- P ) (^ 

^in' 'ppm 


-)^ + p 


I 


IN 


For each assumed pump pressure, the pump flow rate is determined from 
an input data table, equation ( 26 ) and equation (27). The pump outlet 
velocity is then determined from equation (3). 


For each pump flow rate, the quality and flow rate out of each consecutive 
node is a function of the current node pressure, and the node pressure 
is a function of the system pressure dron at the current flow conditions. 
Therefore, for each assumed pump flow condition, an iterative solution 
is required for the system pressure distribution. For the first assumed 
node pressures, P»|^, each node pressure is assumed to be the last 
assumed value for the previous time step. 

A "DO LOOP" is used for consecutive nodes from the pump to the system 
exit to determine the energy and mass balance for each consecutive node. 
Therefore, when entering the "DO LOOP" for a given node, the flow conditions 
and enthalpy leaving the upstream node are known. If the flow ar^a changes 
at the node entrance, the velocity entering the node is determined from 
equation (7). The current vapor volume fraction for the upstream node is 
determined from equation (14) as a function of the current quality and vapor 
density where the vapor density is determined from a thermodynamic table 
as a function of pressure and temperature. The node friction pressure 
drop is then determined from the equations given in paragraph 2.2.1 as a 
function of entering flow conditions. If the flow area changes at the 
node entrance, the entrance pressure loss is determined from equation (18) 
or (19) ard equation (8). 

For the energy and mass balance, the final state is first assumed to be 
two-phase (Xp|\| between 0 and 1) at the current assumed node pressure. 

The fluid temperature is equal to saturation temperature and the wall 
temperature is determined from equation (9). Equations (1), (3), (28), 

(30), (5), (32), and (2) are then solved simultaneously in terms 
of the final mass. This equation is a cubic polynomial. The first 
real, positive root is selected for the final mass, and if no real, 
positive root exists, a real, negative root is selected. After the final 
mass has been determined, the final internal energy, and quality are deter- 
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3.2 (Continued) 

mined from equations (2), (30), and (28). If the final quality from 
the two-phase solution is less than 0, it is set equal to 0, and if the 
final quality is greater than 1, it is set equal to 1. The final density 
IS then determined from equation (30). 


for^finll mL? I two-phase solution is 0 or 1, the above solution 

IZ f internal energy is invalid. If the final quality is 0, 

the final one-phase density is the liquid density, or if the 
IS I , the final one-phase density is the vapor density. The 
density is approximated here from the vapor density table as 


the current assumed pressure and initial temperature, 
final mass is then determined from equation (2). 


final quality 
final vapor 
a function of 


The one-phase 


After the final density has been determined for either the two-phase or 
rat^^exft'^uefni-lfv «lues of density, node exit flow 

(3).’and (2sf ?espectileir determined from equations (4), ( 1 ). 


For the one-phase case (Xp|\j = 0 or 1 ) the final internal energy, effective 
average enthalpy, final temperature, and final wall temperature are 
determined from a simultaneous solution of equations (5), (32), (31) and (8). 

For both the 2-phase and the 1-phase solutions the effective average 
enthalpy is determined from equation (32). 


After the energy and mass balance has been completed for all the system 
nodes at the assumed current pressures, the current system pressure 
distribution is determined from equations (5), (24) and (25)„ Each 
calculated node pressure is compared to the assumed value, and if 
necessary, a new system pressure profile is assumed,. The energy and 
mass balance is repeated at the new assumed pressures until the resulting 
calculated pressures approach the assumed values. 

After the linenode pressure iteration is complete, the calculated pump 
back pressure is compared to the assumed value. If necessary, a new pump 
exit pressure is assumed and the above system energy and mass balance 
and system pressure iteration is repeated at the new pump flow rate. 

After the calculated pump pressure converges with the assumed values, the 
calculations for the time step are concluded. 

3.3 ALTERNATE BRANCH FLOW 


The program has a provision to switch branches of flow to the system 
exit at any prescribed time, e.g. close the vent line and open the 
receiver tank fill line. 
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TRANSIENT CRYOGEN TRANSFER SYSTEM ANALYSIS 

4.0 INPUT DATA 

The Transient Cryogen Transfer Program (TCTP) requires three sets of 
card input data. The first two sets of data are in NAMELIST format 
and the last data set is fixed format input, TCTP also requires a 
data tape (Internal 1/0 Unit 10) and scratch disk space (Internal 
1/0 Unit 11) for cases where the restart capability are exercised. 

4.1 RESTART DATA TAPE INITIALIZATION 

The cryogen transient restart data tape is initialized by writing 
a single 80 byte record containing the integer value - 1000 as the 
first word (4 bytes) in the record. Each of the remaining 19 words 
(76 bytes) should be initialized to zero (or any dummy integer value). 

4.2 NAMELIST INPUT DATA 

The first data set required for TCTP input is the NAMELIST RESTl . 

RESTl contains the indicators, flags and data required to generate a 
restart record on the data tape, and/or to restart a case from the 
data recorded on tape (Restart capability is discussed in Reference 9). 

Variables contained in NAMELIST RESTl are: 

1. RSTFLG - Logical variable used to indicate whether a restart run 

is desired. .FALSE. - Nominal run, no restart from 
existing data. .TRUE. - This run is to be started from 
data stored on tape. 

2. WRSTD - Logical variable used to indicate whether or not a 

restart data file should be written in this run. 

.FALSE, - No data is to be written into the library 
on tape. .TRUE. - Write data into the library on tape. 

3. KEYST - Identification integer for data to be used for a re- 

start. Any positive integer is allowed but it is advised 
that the word be built as follows: 

KEYST = NMMIIJJ 

where N = Record Number for the data (N = 1-9) 

MM = Month (MM = 01-12) 

II = Day (II = 01-31) 

JJ = Year (JJ = 74) 

Thus the third run of June 10, 1974 would have an 
identification integer 

KEYST = 3061074 

NOTE: KEYST card not necessary if RSTFLG = FALSE 
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4. KEYWR 


5. TIMWR 

6. FTIME 

7. VTIME 


Identification integer for data to be written into the 
library on tape. Any positive integer is allowed but 
it is advised that the' word be built as shown in the 
definition of KEYST. 

Time to write the data file for this case into the 
library on tape (SEC), 

Stop time for the current case (SEC). If starting from 
cards, FTIME is overridden by tp in fixed format data. 

Time to switch from vent lines to fill lines in the 
current case (SEC), 


The second data set required by TCTP is NAMELIST INPUT. INPUT con- 
tains variables required to simulate modulating values, pump inlet 
pressures, and storage tank mass levels. The variables in NAMELIST 
INPUT are: 


PUMP/INLET PRESSURE 

APIN(I)* - Array of pump inlet pressures for times given in 

ATPIN, PSIA (maximum of 20 points). 

ATPIN(I)* - Array of times which correspond to the pressures 

in APIN, SEC (maximum of 20' points ) . 

STORAGE TANK SIMULATION 


AZU(I) 

AVZU(I) 

APU(I) 

ATPU(I) 

ATBULK(I) 


- Array of heights in the storage tank referenced to 
2=0. AZU tables points correspond to tank volumes 
in AVZU, FT (maximum of 20 points). 

- Array of volumes in the storage tank. AVZU table 
points correspond to tank heights in AZU, FT**3 
(maximum of 20 points). 

- Array of tank ullage pressures. APU table points 
correspond to times given in ATPU, PSIA (maximum of 
20 points) . 

- Array of times which correspond to ullage pressures 
given in APU, SEC (maximum of 20 points). 

- Array of average bulk temperatures in the storage 
tank. ATBULK table points correspond to times given 
in ATTBLK, DEG R (maximum of 20 points). 


*TCTP gives the user an option of either direct input of pump inlet pressure 

or simulation of the storage tank and feedlines. If data for both simulations 
are included in the NAMELIST data, TCTP automatically defaults to the input 
values of pump inlet pressure versus time. 
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4.2 (Continued) 

ATTBLK - Array of times which correspond to bulk temperatures 

given in ATBULK, SEC (maximum of 20 points). 

NPVZ - Number of table points in AZU. 

ZUI - Initial liquid height in the storage tank referenced 

to z=0 at time zero in the simulation, FT. 

ZTB - Reference elevation of the storage tank exit, FT. 

MODULATING VALVE SIMULATION 


NVALVE(I) 

AVP|3S(J,I) 


ATVP0S(J,I) 

APL0DE(J,I) 

AL0DE(J,I) 


- Array of external node numbers in which modulating 
valves are located (maximum of 5 valves per system 
simulation) . 

- Array of valve positions for each valve corresponding 
to the times given in ATVP0S. J specifies the table 
point, I specifies the valve such that NVALVE(I) is the 
external node number for that valve whose position versus 
time is given in AVP0S(J,I). The range of AVP0S must 
correspond to positions given APL0DE. (Maximum J=20, 
maximum 1=5). 

- Array of times corresponding to the valve positions 
given in AVP0S (maximum J=20, maximum 1=5). 

- Array of valve positions which correspond to the (L/D)p 
given in AL0DE (maximum J=20, maximum 1=5). 

- Array of equivalent length to diameter ratio for the 
valve positions given in APL0DE (maximum J=20, maximum 


WINDMILLING PUMP AND SUCTION LINE SIMULATION 


APMDP(I) - Array of differential pressures (PquT'Pin) Pump flowrate 
calculation during windmilling operation. Pressure data 
must be input in ascending order (PSID) (Maximum of 20 
points) 

APMMDF(I)- Array of pump flowrates corresponding to the differential 
pressures given in APMDP (LBM/SEC) (Maximum of 20 points) 

NPPMDF - Number of table points in APMDP and APMMDF. 

DPMPIN - Diameter of the pump inlet node (FT). 
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4.2 (Continued) 


LPMPIN - Length of the pump inlet node (FT). 

ZPMPIN - Elevation of the pump inlet (FT). 

FPMPIN - E(L/0)g for pump inlet node (actual length/diameter + 

component (L/0)e but not including programmable valve. 

AVPFIN(I)- Array of valve positions corresponding to the times 

given in ATVPIN for valve in the pump inlet node. The 
range of AVPPIN must correspond in magnitude to positions 
given in APLOPN (Maximum of 20 points). 

ATVPIN(I)- Array of times which correspond to valve positions given 
in AVPPIN (SEC) (Maximum of 20 points). 

APLDPN(I)- Array of valve positions which correspond to the (L/D)g 
given in ALDPIN (Maximum of 20 points). 

ALDPIN(I)- Array of equivalent length to diameter ratio for the 

valve positions given in APLOPN (Maximum of 20 points). 

PUMP BY-PASS OR DRAIN LINE FLOW SIMULATION 


BPLODE - E4fo(L/D)e, friction loss coefficient, Kg, for hardware 
elements and actual line length in the BY-PASS line. 

BPA - Area of the BY-PASS line (FT^). 

BPDZ - Differential height (Z^xj-p-Zji^) for the BY-PASS line (FT). 

BPPK - BY-PASS line exit static pressure (PSIA). 

SPECIAL THROTTLABLE BY-PASS LINE SIMULATION 

Special equations were inserted in the TCTP to simulate throttlable 
branch flow from a system node to the supply tank reservoir (KSC LOX 
by-pass) as a function of the system node exit pressure. These equations 
are labled "KSC BY-PASS FLOW" in the program, subroutine CALQL8, 
listing with comment cards. If this simulation is not used, these 
equations should be removed from the program, or the input data should 
be set such that the calculated by-pass flow is zero. The required 
NAMELIST data input for this simulation are: 
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4.2 (Continued) 

ACBPl - Array of by-pass line flow coefficients (288g A^p/Kr) 
corresponding to the times given in ATCBPl ^ ^ ^ 

(Maximum of 20 points). 

ATCBPl - Array of times corresponding to the by-pass line flow 
coefficients (Maximum of 20 points). 

NPCBPl - Number of table points in ACBPl. 

Ifthis simulation is used, the internal node numbers at the by-pass 
exit and the next two downstream nodes and the differential height 
(ZeXIT-Zin) for the by-pass line must be inserted directly in Subroutine 
CALQL8 code to calculate the by-pass flow as discussed in paragraph 

2.4.2. 

SPECIAL SYSTEM BLEED FLOW SIMULATION 

Special equations were inserted in the TCTP to simulate a bleed flow 
from the system at a given node exit as a function of the mode exit 
pressure. These equations are labeled "KSC BLEED FLOW", "SPECIAL 
EQUATIONS TO SIMULATE ENGINE BLEED", "BLEED FLOW", and "KSC BLEED 
FLOW" in the program listing (Subroutine CALQL8). If this simulation 
is not used, these cards should be removed from the program, or by- 
passed with bleed flow set to zero. If this simulation is used, the 
following must be inserted directly into Subroutine CALQL8 code to 
calculate the bleed flow as discussed in paragraph 2.4.2. 

(1) Internal node numbers at the bleed flow exit and the next 
two downstream nodes. 

(2) Bleed flow exit (vent line) pressure. 

(3) Bleed flow system differential height constant (Z£xjj-Zjj^)/144. 

(4) Bleed flow system flow coefficient (288 g^ ^bL'^'^E^’ 

PUMP PERFORMANCE SIMULATION OPTION 

An option has been incorporated in the program to compute required 
pump speed as a function of input system entrance flow rate and 
computed pump discharge pressure instead of computing system entrance 
flow rate as a function of input pump speed and computed pump discharge 
pressure (paragraph 2.4.1). In order to excercise the option to input 
pump flow rate instead of pump speed, the required NAMELIST data are: 
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4.2 (Continued) 

MDPUMP = 1 Option to input pump discharge flow rate as a 
function of time. The program will determine 
required pump speed as a function of pump flow 
rate and pressure rise, 

MDPUMP = 0 Option to input pump speed (RPM) as a function 
of time. The program will determine pump flow 
rate as a function of pump speed and pressure 
rise. 


AMDFI Array of pump discharge flow rates corresponding 

to the times given in ATMDFl (Maximum of 20 points). 

ATMDFl Array of times corresponding to the pump flow 

rates (maximum of 20 points). 

If no values of pump flow rate are input, the program will default 
to MDPUMP = 0. 


4.3 FIXED FORMAT INPUT DATA 

The fixed format input data are given in three logical sections and 
set UD in the following order: 

Section 1 System Constants 

Section 2 Nodal Data 

Section 3 Tables 


Unless otherwise noted, all values are floating point values and may be 
keypunched anywhere in the indicated card field (columns). 

Integer values must be right adjusted in their fields. 

A blank field (no punches in the field) is equivalent to a zero (inteoer 
or floating point). ^ 

4.3.1 Section 1 - System Constants 

CARD COLUMNS CONTENTS 

^ -one to seventy alphanumeric characters 

(including blanks) to be used as a print title 
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CARD 

2 


3 


4 


! 


1 


T 

! 




(Continued) 

COLUMNS 

1-10 

11-20 

21-30 

31-40 

41-50 


1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

1-10 

11-20 

21-30 

31-40 

41-50 

51-60 


1-10 


CONTENTS 

At - time increment (sec) 

tp - final time (sec). If restarting from tape, 
tp is overridden by FTIME NAMELIST data. 

tv - vent close time (sec) 

Atp - print time increment (sec) 

print option (integer) 

= 0 - do not print optional data 
= 1 - print optional data 

BLANK SPACE - NOT USED 

g - system acceleration (ft/sec2) 

yy - vapor viscosity (lb(mass)/ft-hr) 

yj_ - liquid viscosity (lb(mass)/ft-hr) 

Ky - vapor thermal conductivity (BTU/ft-hr-®R) 

K|_ - liquid thermal conductivity (BTU/ft-hr-°R) 

~ liquid bulk modulus (psi) 

Ppy - vapor Prandtle number 
PrL - liquid Prandtle number 

CpL - liquid specific heat (BTU/LB(Mass)“R) 

Cpy - vapor specific heat at constant pressure 
(BTU/Lb(Mass)°R) 

Cyw - vapor specific heat at constant volume 
(BTU/lb(mass)°R) 

Np - pump external node number (integer) 


11-20 

21-30 

31-40 

41-50 

51-60 

61-70 


Nj - transfer line low and high external node 
numbers (integers) 

Nw - vent line low and high external node 
numbers (inteaers) 

Np - fill line low and high external node 
numbers (integers) 
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4.3.1 (Continued) 

CARD COLUMNS 


6 1-10 
11-20 
21-30 
31-40 
41-50 
51-60 
61-70 


CONTENTS 


external node numbers (integers) of the nodes 
for which print is desired. One to seven node 
numbers per card (zero node numbers are ignored). 
Repeat this card until all nodes are given. 

To print all nodes omit this card(s). (Maximum 
of 50 input values). 


7 


4.3.2 


1-70 


Blank - a blank card to signal the end of print 
nodes input. 


Section 2 - Nodal Data 


CARD COLUMNS 

1 1-10 

11-20 
21-30 
31-40 
41-50 

51-60 

61-70 


2 1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

61-70 


3 1-70 


is 

QtLlLiTy 


CONTENTS 

N - external node number (integer) 

D|^ - diameter (ft) 

L|^ - length (ft) 

Z|^ - elevation (ft) 

Pext " external heat rate (BTU/ft^sec) at Tj^j-T.| 

- Loss coefficient for area change 

T|^j - initial temperature (°R) if a zero or 
negative value is given, the program 
will calculate saturated Twt as a function 
of pr,,„ using Table 2 

A|M|, - Actual ''nside heat transfer surface area 
if different from 

P|\|j - initial pressure (psi(abs)) 

Xf^j - initial bulk quality 

z(L/D) for node(actual length/diameter + component (L/D) 
^W|\| * * wall specific heat 

MDFjyj - initial flow rate, Input -1 if program 
calculation is desired. 

TWj^jj - Initial wall temperature. 

Input (-1) if program calculation is 
desired for steady state initial 
conditions. 

blank - a blank card to signal the end of 
nodal data 
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4.3.2 (Continued) 

fast^nodr^ ^ i^epeated for each node. Card 3 follows card 2 of the 

not have to be given in any order. The program will sort 
them byline section, and by node number, and assign internal node 
sequencing. 

No external node number may be zero or negative. 

4.3.3 Section 3 - Tables 

The following tables are used in the TRANSIENT CRYOGEN TRANSFER PROGRAM. 
Table #1 Vapor density as a function of pressure and temperature 
PV = (P> T) 

Table §2 Saturation temperature as a function of pressure 
Tj = fg <P> 

Table #3 Vapor saturation internal energy as a function of temperature 

%s = fa <fs> 

Table #4 Saturated liquid internal energy as a function of temperature 
Ufs " '4 (f) 

Table #5 Friction factor correction as a function of vapor volume 

<f?) = fg ('-«) 

Table #6 Single phase friction factor as a function of two phase 
Reynolds' number 

(fOp) = fg (^) 

Table #8 Exit node pressure as a function of time 

Pk ' fs (f' 
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4.3.3 (Continued) 


Table #9 (Pump flow rate/RPM) as a function of (pressure/’RPM^) 


Table #10 


( 


JL) 

RPM^ 


Rpm > 


Pump RPM as a function of time 
RPM = f^Q (t) 


Table #11 


Pump discharge temperature as a function of time - Input 1 
if pump is "windmilling" 



(t) 


The table input cards have the following columns and content. 
CARD COLUMNS DATA 

3.1 1-10 Table numbers 


n-20 Tg - Built in table type - Tables #1 to #4 

have been built into the program for 
various types of fluid. Tables #5 and #6 
have been built in also. When a built in 
table is to be used, Tr is a number defininq 
the types of fluid. (Tables #1 to 4) 

Tg = 1 - LOX 

Tg = 2 - Liquid H 2 (future addition) 

For tables #5 and #6 Tg = 1 

NOTE: set Tg = 0 to read the table from cards. 

21-30 - BLANK SPACE - NOT USED 


31-40 N" - blank space - NOT USED 


41-50 


X values scale factor 
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(Continued) 

CARD COLUMNS DATA 


51-60 Y values scale factor 

61-70 Z values scale factor 

When the table is given on cards, it may be convenient to scale the values 
for two dimensional tables, the X scale factor scales the independent 
values and the Y scale factor the dependant variables. For three dimensional 
the X and Y scale factors scale the independent variables, the Z scale factor 
the dependent variable. In all cases the scaling is (desired value) = 

(scale factor) times (input value). For tables #1 to #6 when Tg 
is given. Card 3.1 is all that is required as input. 


3.2 


For two dimensional tables only (tables #2, #4, 
#5, #6, #7 (all nodes), #8, #9, #10, and #11) 


11-20 

21-30 

31-40 

41-50 

51-60 

61-70 


3.3 to 
3.K-1 


X-j - the independent value 
Y-j - the dependent value at Xi 


NOTE: The values are given in pairs. If the 

card columns for any pair are left blank, or 
both values are zero, then the pair is dis- 
carded. To obtain a true zero-zero pair, 
input one value very small (i.e. l.OE - 60). 
At least one paiv must be given on each card. 


Repeat card 3.2 as necessary to input the full 
table 


I 
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4.3.3 (Continued) 


CARD 

COLUMNS 

DATA 

3.K 


One blank card 

3.2 


For three dimensional tables only (tables #1 
and #3) 


11-20 

X-| - first independent value 


21-30 

Y-| - second independent value 


31-40 

- dependent values at X-j , 


41-50 



51-60 

h 


61-70 

h 



NOTE: The values are qiven in sets of three. 

If the card columns for any set are left blank 
or all values of the set are zero then the set 
is discarded. To obtain a true zero-zero-zero 
set, input one value very small (i.e., 1.0 E-60) 
At least one set must be qiven on each card. 

3.3 

to 

K 


Repeat card 3.2 as necessary to input the full 
table 

3K 


One blank card 

Card 3.1 must be given 
#3, #4, #5, or #6 the 
fluid other than LOX, 

for each table. If not given for tables #1, #2, 
built-in table(s) will be used. To default to a 
use a card 3.1 with the table numbers = 0, and 


the types of fluid (Tb) set for the fluid desired. Table #7 is not 
needed for all nodes. For the nodes for which a table #7 is not given 
a value of Kg = 0 will be used at computation time. 


Cards 3.2 to 3.K must be given for each table unless the value of Tb 
( tables #1 to 6 and #12) or the value of N" (table #7 only) indicate 
otherwise in which case cards 3.2 to 3.K must be omitted. 

Card 3.K is required for all tables when the values are given as input 
(card 3.2 etc.). Tables may be given in any order. 
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4.3.4 Deck Se^iiiencing 

Card Columns 71 to 80 are not used by the program, and therefore can be 
used for deck sequencing. A suggested method of sequencing is as 


follows; 



COLUMNS 

VALUE 


72 

Section number 


74-75 

Table number (section 

3 only) 

77-79 

Node numbers (section 

2) 

80 

Card number 


4.4 

PARALLEL AND BRANCH FLOW VERSION OF TCTP 


During the SSME Feasibility Study (Reference 11), a special version of 
TCTP was developed to simulate systems with multiple flow paths. This 
version of TCTP will incorporate either oxygen or hydrogen as the 
working fluid by setting the fluid type indicator to one or two, 
respectively, and automatically processes multiple flow paths if the 
input data are arranged according to the rules set forth in this 
section. 


The following definitions are basic to the multiple flow path version 
of TCTP: 


Main Flow Leg 


Parallel Flow Leg 


Branch Flow Leg 


An arbitrarily chosen flow path through a 
system. This path must start at the system 
supply point (usually a pump) and terminate 
at a plenum which is outside the system. 
The pressure in the plenum may be constant 
or time variant. 

A flow path which emanates from a point in 
the main flow leg and terminates into a 
different point in the main flow leg. 

A flow path which emanates from a point in 
the main flow leg and terminates into a 
plenum. The pressure in the plenum may be 
constant or time variant, but the plenum 
must be external to the system. 
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Junction points of the main flow leg and branch 
leg entrances and/or parallel leg entrances 
and exits. Mixing regions are always located 
at the downstream (away from the system supply 
node) end of a main flow leg node. One mixing 
region may contain one or a combination of 
entrances and/or exits. 

The following rules must be adhered to in formulation of the system 
description for input to TCTP: 

1) The maximum number of legs' will be six per system: one main flow 

leg and any combination of five (or less) parallel and branch 
legs. The minimum system description shall contain at least one 
leg: the main flow leg. 

2) All branch and parallel flow legs will emanate from the main flow 
leg; parallel legs will terminate into the main leg. 

3) No nested parallel flow paths are allowed; each parallel leg must 
be terminated before another parallel leg is initiated along the 
main leg. 

4) Branch and parallel legs may be initiated at any node in the main 
flow leg with the exception of the system supply node and the system 
exit node (first and last nodes). 

5) A minimum of two nodes shall be placed in each parallel or branch 
leg. 

The nodal data required for the physical description of each node in 
the system are identical to those given in Section 2 (cards 2.1 and 
2.2). However special instructions are required in assigning node 
numbers for the system. The node numbers must be assigned in 
ascending order from the system supply node to the main leg exit, 
then node numbers for the first leg in ascending order, node numbers 
for the second leg, etc. All numerical values assigned as node numbers 
are arbitrary but must be assigned in ascending values. Leg node 
numbers must always be larger than the main leg exit node numbers 
and must be grouped by legs. Numerical values for leg node numbers 
must fall within the "vent" range as identified by card 1.5. 

The following variables are input in NAMELIST format using NAMELIST 
INPUT: 

NLEGS Number of flow legs in the system including the main flow 
leg. (NLEGS. GE.l) 


S' 

i 


4.4 (Continued) 
Mixing Region 
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4.4 (Continued) 


LENT(I) 

LEXT(I) 

LEG(I) 

ATPK(J,I) 

APK(J,I) 

NPPK(I) 


Node number for the main flow leg node which is the 
supply node for a parallel or branch leg. (NOTE: 
legs always exit from .the downstream end of a main 
flow leg node). LENT(l)=the pump node number, I=leg 
number. 

Node number for the main flow leg. node into which a 
parallel leg terminates. LEXT(l)=the last node in the 
main flow leg. LEXT(I) is input .LE.O for branch flow 
legs. I=leg number. 

Node number of the first node in each leg of a system. 
LEG(l) = pump node number. LEG(I), I.GT.l, should 
always be greater than LEXT(l) and all node numbers 
for leg I should be less than LEG(I+1); node numbers 
in the last leg node should be less than NVENT(2) 

(card 1.5). All LEG(I) values must fall into the 
range NVENT(1 )<LEG(I)<NVENT(2) for I>1 . I=leg number. 

Array of times corresponding to pressures given in 
APK(J,I) (seconds). I=leg number, J=table point 
number (maximum of 20 points per leg). 

Array of plenum pressures for the main leg and all 
branch legs corresponding to the times given in ATPK 
(J,I). (PSIA) I=leg number, J=table point number 
(maximum of 20 points per leg). 

Number of points in the APK(J,I) array. 1=1 eg number. 


Mixing regions are automatically assigned for each leg entrance and 
exit. (No mixing region is assigned for leg one). Three coordinates 
are defined at each mixing region as follows: 


SYSTEM 

SUPPLY 



SYSTEM 

EXIT 
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4.4 (Continued) 

These coordinates are used to identify flowrate, enthalpy and 
density crossing the boundaries of the mixing region at a leg 
entrance or exit. When combinations of entrances and/or exits 
exist at mixing region, initial values for flowrate at boundaries 
(1) and (3) are input as total values for each leg. The following 
variables are input via NAMELIST INPUT and conform to this coordinate 
convention: 


MDFSP(K,I) 

HPTSP(K,D* 

RH0SP(K,I) 

MDFSPI(K,I) 

HPTSPKK,!)* 

RH0SPI(K,I) 


- Initial flowrates crossing boundaries of a mixing 
region for a parallel or branch leg entrance. Flows 
are positive left to right for main leg nodes. Flows 
are positive entering a leg; negative if the> are 
being expelled from the leg at the entrance. (LBj^/ 
SEC). I=leg number, K=coordinate. 

- Initial Enthalpy of fluid crossing boundaries of a 
mixing region for a parallel or branch leg entrance 
(BTU/LB^^). I=leg number, K=coordinate. 

- Initial density of fluid crossing boundaries of a 

mixing region for a parallel or branch leg entrance 
(LB|^/FT3). 1=1 eg number, K=coordinate. 

- Initial flowrates crossing the boundaries of a 
mixing region for a parallel leg exit. Main leg 
node flows are positive left to right. Flow being 
expelled from a parallel leg into the mixing region 
is positive. (LBn^/SEC). I=leg number, K=coordinate. 

- Enthalpy of fluid crossing boundaries of a mixing 
region for a parallel leg exit (BTU/LBn^). I=leg 
number, K=coordinate. 

- Density of fluid crossing boiindaries of a mixing 
region for a parallel leg exit (LB|^/FT3). I=leg 
number, K=coordinate. 


♦Enthalpy for hydrogen mi’st be defined in a manner consistent with 
NBS-TN-617 (Reference 12' . Oxygen enthalpy must be defined on a 
base consistent with NASA SP-3071 (Reference 13). 
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4.4 (Continued) 

An example of a system with parallel and branch legs is shown 
schematically in Figure 6. A card image of the input data deck 
is given in Figure 7. 
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F.P.B. 

MANIFOLD 



PK(1) = OXYGEN PREBURNER MANIFOLD 
PRESSURE 

PK(2) LPFT INLET PRESSURE 
PK(3) = 0.0 (PARALLEL LEG) 

PK(4) = FUEL PREBURNER MANIFOLD PRESSURE 


KEY: 



NODE 




VALVE 


FIGURE 6 


SSME HYDROGEN SYSTEM SCHEMATIC WITH TCTP NODES IDENTIFIED 
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FIGURE 7 CARD IMAGE OF INPUT DATA FOR SSME HYDROGEN SYSTEM 
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FIGURE 7 CARD IMAGE OF INPUT DATA FOR SSME HYDROGEN SYSTEM (CONTINUED) 
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ATTACHMENT 2 


STEADY STATE CRYOGEN TRANSFER SYSTEM ANALYSIS 


1.0 CONTROLLING EQUATIONS 


Many of the equations used in the steady state program are identical to 
those given in Attachment 1 for the transient program. Therefore, for 
convenience the additional steady state equations will be numbered 
consecutively with the equations in Attachment 1. 


For steady state conservation of energy, equation (4) reduces to. 


h- = h- + ^'n-1 , Vn-1 

^ N ^ N-1 ^ gcU 2gcJ 


•^EXT 




V|\| 




(33) 


For conservation of momentum, equation (5) reduces to, 

.-2 .... 2 

TWgT " PN^^Tzfe ^^N " Vl^^ " ^''fN. 


^N “ '"n-1 


^N-1 ''N-1 " Pn'''^ 
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STEADY STATE CRYOGEN TRANSFER SYSTEM ANALYSIS 
2.0 SUBSIDIARY EQUATIONS 

In addition to the fluid property equations given in Attachment 1, the 
one-phase (X[^ = 0 or 1) node temperature is determined as a function 
of enthalpy, 


T _ -r '^S " ^N 

'n 'sn ■ C 

(35) 

where h^ = h^^ and C = Cpy 

for X^ = 1 

and h^ = h|_g and C = Cpj^ 

for X,^ = 0 

The node quality is determined as a 
assumed node pressure. If the node 
liquid and vapor values (hj^g ^ 

function of the exit enthalpy at the 
exit enthalpy is between the saturated 
hys) 

V _ ^N ■ ^LS 

^ ^vs " \s 

(36) 

^ *^LS’ ^N ° 

(37) 

and if hj^j > h^^, = 1 .0 . 

00 

oo 

• 
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STEADY STATE CRYOGEN TRANSFER SYSTEM ANALYSIS 


3.0 SYSTEM ANALYSIS 

After the known system conditions and geometry are established, the 
procedure outlined on Figure 8 is utilized to determine the steady state 
!v w rate corresponding pump RPM, node properties and pressure gradient. 

L" finally, the program computed both the node bulk average quality, 

Xm, and the local quality leaving the node, Xj^, and corresponding 
vl!iues of density and vapor volume fraction. However, in the current 
program, the node bulk and exit qualities are considered to be equal. 
Therefore, the redundant equations for bulk and local values of density 
and vapor volume fraction could be eliminated with some program nomen- 
clature changes. 

The pump discharge flow entering the transfer system is assumed to be 
subcooled at the input values of temperature and pressure. 

The pump exit density is set equal to the saturated liquid density for 
the given cryogen, and the saturated internal energy is determined from 
an input table as a function of input temperature. The pump exit 
enthalpy is determined from equation (32). The quality and vapor volume 
fraction entering the transfer system are set equal to zero, and the 
velocity is determined from equation (3) for the assumed flow rate. 

Since the hydrostatic head, two-phase frictional pressure drop, and 
momentum flux are functions of the fluid quality and the quality is a 
function of flow rate, an iterative solution is required to determine 
the steady state flow rate for the known entering and exit node pressures 
and system external heat flux. 

A "DO LOOP" is used for consecutive nodes, 2 through K where K is the 
system exit node, to determine each node pressure, thermodynamic state, 
and exit velocity. Therefore, when entering the "DO LOOP" for a given 
node, the flow conditions and enthalpy leaving the upstream node are 
known. If the flow area changes at the node entrance, the velocity 
entering the node is determined from equation (6) and the entrance 
pressure loss (Pr, i - Pm i) is determined from equation (18) or (19) 
and equation (s)* The vapor volume fraction for the upstream 

node is determined from equation (14) as a function of the quality and 
vapor density where the vapor density is determined from an input table 
as a function of pressure and temperature. The node friction pressure 
drop is then determined as a function of entering flow conditions from 
the equations given in paragraph 2.2.1, Attachment 1. 

Since the energy balance and pressure gradient across each node is a 
function of the node pressure and density, an iterative solution is used 
to determine the pressure and density. For the first assumption, the node 
density and pressure are set equal to the upstream node values. The 
saturated values of temperature, liquid internal energy, and vapor 
enthalpy are determined from input thermodynamic tabies at the assumed 
node pressure. The saturated liquid enthalpy is determined from 
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3.0 (Continued) 

equation (32) and the node exit velocity is determined from equation (3) 
for the assumed flow rate and node density. The node exit enthalpy 
is then determined from equation (33). 


The node quality and temperature are determined as a function of the node 
exit enthalpy at the assumed pressure. The quality is determined from 
equation (36), (37) or (38). If the node exit enthalpy is between the 
saturated liquid and vapor values, the node temperature is equal to 
the saturation value at the assumed pressure. If the node is super- 
heated or subcooled, the temperature is determined from equation (35), 

The vapor density is determined from an input table as a function of 
temperature at the assumed pressure, and the node density is determined 
from equation (30). If the assumed density is not within 5 percent of 
this calculated value, a new value is assumed and the above node velocity 
and energy calculations are repeated. If the assumed and calculated 
densities converge, the velocity is corrected to the last calculated 
density, and the pressure change across the constant area node is determined 
from equation (34). The node exit pressure is then determined from. 


+ (P^ 


- P' 


) + (P 


N-l 


- P 


N-l) 


If the calculated node pressure is not within 0.01 psi of the assumed value, 
a new pressure is assumed and the above energy/momentum balance is 
repeated. After the "DO LOOP" is completed through the system exit 
node, the calculated system exit node pressure is compared to the known 
exit pressure. If the calculated system exit pressure is not within 
0.01 psi of the known system exit pressure a new flow rate is assumed. 

The above pump velocity calculation and "DO LOOP" is repeated for the 
entire system with new flow rate assumptions until the calculated 
system exit pressure converges with the known value. 

After the system flow rate is determined for the known system entrance 
and exit pressures, the two-phase choked flow relation derived in 
paragraph 4.2, Attachment 1, is used to determine the choked flow rate 
for each node. If the previously calculated flow rate exceeds the 
choked flow rate for any node, the node "NO" with the highest positive 
value of (I^ - 1^*) is identified. The system flow rate and pressures 
are then adjusted according to the procedure outlined in Figures. 

For the selected example problem, choked flow was not encountered. 

Therefore, the choked flow portion of the program has not been checked 
out. 


After the steady state flow rate is determined, an iterative procedure 
is used to determine the corresponding pump RPM at the known pump dis- 
charge pressure. This RPM is a required input for the transient program 
to determine the initial flow rate. The first assumption for RPM must 
be input to the steady state program. For each assumed RPM, the pump 
flow rate, t^r, is computed from the input pump performance table 
(equation(26) and equation (27)). If the computed pump flow rate is not 
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3.0 (Continued) 

within 0.01 LBS/SEC of the system f1o.w rate, a new RPM value is assumed 
and the pump flow rate calculation is repeated. After the pump RPM 
corresponding to the steady state system flow rate is determined the 
steady state program calculations are terminated. 
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FIGURE 8 STEADY STATE CRYOGEN TRANSFER PROGRAM LOGIC 
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FIGURE 8 STEADY STATE CRYOGEN TRANSFER PROGRAM LOGIC (Continued) 
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STEADY STATE CRYOGEN TRANSFER SYSTEM ANALYSIS 

4.0 INPUT DATA 

The Steady State Cryogen Transfer Program (SSCTPl requires two sets of 
card input data. The first data set is in NAMELIST format and the second 
set is fixed format input. SSCTP requires no tape or disk allocations. 

4.1 NAMELIST INPUT DATA 

The first data set required for SSCTP input is NAMELIST WALC0N. WALC0N 
contains data which is punched out for use in the TCTP. The variables 
contained in NAMELIST WALC0N are: 

WALLM(N) - node wall mass (LB(MASS)) 

WALLH(N) - node wall specific heat (BTU/LB(MASS)°R) 

MDF(N) - mass flowrate out of node N (LB(MASS)/SEC) 

P(N) - pressure at exit of node N (PSIA) 

X(N) - bulk quality in node N 

XP(N) - local quality in node N 

T(N) - temperature in node N (°R) 

4.2 FIXED FORMAT INPUT DATA 

The fixed format input data are given in three logical sections and set 
up in the following order. 

Section 1 - System Constants 
Section 2 - Nodal Data 
Section 3 - Tables 

Unless otherwise noted, all values are floating point values and may be 
key-punched anywhere in the indicated card field (columns). Integer 
values must be right adjusted in their fields. A blank field (no punches 
in the field) is equivalent to a zero (integer or floating point). 
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4.2.1 

Systan Constants 

CARD 

COLUMNS 

CONTENTS 

1 

1-70 

title to be used as print title 

2 

1-10 

print option (Integer) 

= 0 do not print optional data 
= 1 print optional data 

3 

1-10 

Pamb - ambient pressure (PSIA) 


11-20 

q - system acceleration (FT/SEC ) 


21-30 

Pl - liquid density (LB(MASS)/FT^) 


31-40 

liy - vapor viscosity (LB(MASS)/FT-HR) 


41-50 

M| - liquid viscosity (LB(MASS)/FT-HR) 

4 

1-10 

- liquid bulk modulus (PSI) 


11-20 

Cpi^- liquid specific heat BTU/LB(MASS)“R) 


21-30 

Cpw- vapor specific heat at constant pressure 
(BTU/LB(MASS)°R) . 


31-40 

C.,.,- vapor specific heat at constant volume 
(BTU/LB(MASS)°R) 

5 

1-10 

_ pump outlet temperature (°R) 


11-21 

P|^j- exit node pressure (PSIA) 


21-30 

M - assumed system flow rate 


31-40 

P(l) - pump discharge pressure (PSIA) 


41-50 

RPM - Pump RPM Assumption (RPM) 

6 

1-10 

Np - pump external node number (integer) 


11-20 

Nj - transfer line low and high 


21-30 

external node numbers (integer) 


31-40 

Ny - vent line low and high 


41-50 

external node numbers (integer) 
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4.2.1 

(Conti nued) 


CARD 

COLUMNS 

CONTENTS 


51-60 

Np - fill line low and hiqh 


61-70 

external node numbers (inteqer) 

7 

1-10 

(Np) - external node numbers (integer) of 


11-20 

the nodes for which print is desired, 


21-30 

One to seven node numbers per card (zero 


31-40 

node numbers are ignored). Repeat the 


41-50 

card until all nodes are given. To print 


51-60 

61-70 

all nodes, omit this card(s). 

8 

1-70 

Blank - a blank card to signal the end of 
print nodes input. 

4.2.2 

Section 2 - 

Nodal Data 

CARD 

COLUMNS 

CONTENTS 

1 

1-10 

N^i - external node number 


11-20 

D|^ - diameter (FT) 


21-30 

L^^ - length (FT) 


31-40 

- elevation (FT) 


41-50 

Pext " external heat flux (BTU/FT^SEC) 


51-60 

- loss coefficient for area change 


61-70 

Option Constant (integer) 

= 0 use to calculate friction pressure 

= 1 use 

2 

11-20 

r(L/D) for node (actual length/diameter + 
component (L/D)e) 

3 

1-70 

Blank - a blank card to signal the end of 


nodal data, 
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f 1 


4.2.2 (Continued) 

Card 1 is repeated for each node. Card 2 follows the definition of the 
last node. 

The nodes do not have to be given in any particular order. The nroqrari 
will sort them by line section and node number and assign internal node 
sequencing. 

No external node number may be zero or negative. 


4.2.3 Section 3 - Tables 

Table #1 Vapor density as a function of pressure and temperature 
pv = f^ (P, T) 

Table #2 Saturation temperature as a function of pressure 
Ts = f2 (P) 


Table #3 


Table #4 


Table #5 


Liquid saturated internal energy as a function of temperature 
(Table #4 in transient program) 


U 


LS 


fo (T) 


Friction factor correction as a function of vapor volume 
(Table #5, transient program) 

"’*^4 0 " “) 

Single phase friction factor as a function of Reynolds number 
(Table #6 in transient program) 

(fo) “ O 


p 

Table #7 (Pump flow rate/RPM) as a function of (Pressure/RPM ) 
(Table #9 in transient program) 

( ^1 


£ /P-Pambx 
T? \ 

' RPM^ 
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4.2.3 (Continued) 

Table m Saturation vapor enthalpy as a function of tomperaturo 
(Table ^]2 in transient proqram) 

hvs = fa <f) 

The table input cards have the followinq columns and content. 


CARD COLUMN CONTENTS 

3.1 1-10 Nj - table number (inteqer) 

n-20 NTi^ - Built in table type (integer) - Tables 

#2, #3 and #8 have been built into the program for 
various types of fluid. Tables m and #5 have 
been built inalso. When a built in table is to 
be used, NTj^ is a number defining the type of 
fluid (Tables A1 , #2, #3, and #8 . 

NT(j = 1 - LOX 

NT(-! = 2 - LH 2 (a future addition) 

For Tables #4 and 45 NTj^ = 1 

NOTE: Set NTp = 0 to read the value from cards. 

% - ^Jsed for input of Table 46 and is the node 
number for this Table 46 (inteqer) 

31-40 N^ - (Used for Table 46 only. Table 46 for node 

Nm (see columns 21-30 above) is identical to 
Table 46 for node N|,j. 

41-50 X values scale factor 

^1“50 Y values scale factor 

61-70 Z values scale factor 

NOTE: Right adjust exponential scale factor. 

When the table is given on cards, it may be convenient to scale the values. 
For two dimensional tables, the X scale factor scales the independent 
variables. For three dimensional tables, the X and Y scale factors scale 
the independent variables, and the Z scale factor the dependent variable. 

In all cases the scaling is (desired value) = (scale factor) times (input 
value). 

NOTE: For Tables #1 - 45 and 48, when NT(j, is given, and for Table 46 

when N,^ is given (non-zero). Card 3.1 is all that is required as input. 


21-30 


31-40 


41-50 

51-60 

61-70 
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4.2.3 (Continued) 

CARD COLUMNS CONTENTS 

3.2 -- For tv/o dimensional tables only (Tables «2, 

n, M, #5, #6 (all nodes), #7 and #8) 

11-20 X-| - the independent value 

21-30 - the dependent values at X^. 

31-40 Xg 

41-50 Y 2 

51-60 X 3 

61-70 Y 3 

NOTE: The values are given in pairs. If the card columns for 

any pair are left blank, or both values are zero, then the 
pair is discarded. To obtain a true zero-zero pair, input 
one value very small (i.e. 1.0 E -6.0). At least one 
pair must be given on each card. 

3.3 - Repeat card 3.2 as necessary to input the 

to complete table. 

3.K-1 


3.K 

3.2 


- One blank card to signal the end of the table. 

For three dimensional tables only (Table #1) 

11-20 X.J - first independent value 

21-30 Y-j - second independent value 

31-40 Z-j - dependent value at X^. , Y^ 

41-50 Yg 

51-60 Y^ 

61-70 Z 2 

NOTE: The values are given in sets of three. If the card columns 

for any set are zero, the set is discarded. To obtain a 
true zero-zero-zero-set, input one value very small 
(l.OE-6.0). At least one set must be given on each card. 




I 
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4.2.3 (Continued) 


CARD COLUMNS C ONTENTS 

Repeat card 3.2 as necessarv to input complete 
table. 


- One blank card to signal the end of the table. 

Card 3.1 must be given for each table. If not given for Tables #1, #2, 

#3, #4, #5, or #8, the built-in table(s) will be used. To default to 
a fluid other than LOX, use a card 3.1 with the table number = 0, 
and the type of fluid (NTg) set for the fluid desired. Table #6 is not 
needed for all nodes. For the nodes for which a table #6 is not given, 
a value of Ke = 0.0 will be used. 

Cards 3.2 to 3.K must be given for each table unless the value of NTb 
( tables #1 to #5 and #8) or the value of N|;l (table #6 only) indicate 
otherwise, in which case cards 3.2 to 3.K must be omitted. 

Card 3.K is required for all tables when the values are given as input 
(card 3.2 etc.). 

Tables may be input in any order. 


3.3 

to 

3.K-1 


4.2.4 Deck Sequencing 

Card columns 71 to 80 are not used by the program, and therefore can 
be used for deck sequencing. A suggested method of sequencing is as 
follows: 


COLUMNS 

72 

73-74 

75-78 

79-80 


VALUE 

Section number 

Table number (Section 3.0 only) 

Node number (Section 2,0 and Section 3.0 
tables #6 only) 

Card number 
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ATTACflf'’f;NT 


1.0 CRYOPEN TRANSFER NOMENCLATURE 
SYMBOL DEFINITION 

A Area 

C Specific Heat 

D Diameter 

Eu Euler Number 

F Force 

Tq Sinqle Phase Friction Factor Evaluated 

at the Mixture Reynolds Number 

fjp Two - Phase Friction Factor 

q Gravitational Acceleration or Orbital 

Transfer System Acceleration 

9|- Constant of Proportional ity 

h Node Bulk Enthalpy 

h" Average Enthalpy at a Node Exit 

h^ Heat Transfer Film Coefficient 

J Conversion Factor 


K Friction Loss Factor 

k Thermal Conductivity 

1 Line Segment(or node) length 

M Mass in the Control Volume or Node 

M Mass Flow Rate 

N Pump Speed 

N^j Nusselt Number 


UNITS 

Ft^ 

BTII/l.B^^"P 

Ft 

Dimensionless 

LBp 

Dimensionless 

Dimensionless 

Ft/Sec^ 

LBp-SEC'^ 

BTlI/LBp, 

BTU/LB^^ 

btli/ft2-hr-°r 

FT-LBp 

"BTU 

Dimensionless 

BTU 

FT-HR-°R 

FT 

Lf^M 

lb^/sec 

RPM 

Dimensionless 
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1.0 

(Continued) 


SYMBOL 

DEFINITION 

UNITS 

P 

Static Pressure 

PSIA 

P' 

Static Pressure Downstream of an 
Area Change 

PSIA 

Pr 

Prandtle Number 

Dimensionless 

q 

Heat Transfer Rate 
(BTU/FT2-SEC in computer program) 

BTU/SEC 

T 

Temperature 

°R 

t 

Time 

SEC 

u 

Internal Energy 

BTU/LB^ 

V 

Volume 

Ft3 

V 

Velocity 

FT/ SEC 

X 

Node Bulk Quality 

Dimensionless 

X" 

Average Quality at a Node Exit 

Dimensionless 

Z 

Elevation at a Node Exit 

FT 

a 

Node Bulk Vapor Volume Fraction 

Dimensionless 

a" 

Average Vapor Volume Fraction at a 
Node Exit 

Dimensionl ess 

VI 

Viscosi ty 

LB^^/FT-SEC 

P 

Node Bulk Density 

lb,^/ft^ 

o' 

Average Density at a Node Exit 

LBf,/FT^ 

^tt 

Marti nelli Two Phase Parameter 

Dimensionless 

SUBSCRIPTS 

DEFINITION 


1 

Pump Discharge 


A 

Area Change 


C 

Calculated 


EXP 

Experimental 


E 

Equivalent Value 
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1.0 (Continued) 
SUBSCRIPTS 


F 

I 

IN 

L 

N 

0 

S 

TP 

V 

W 


DEFINITION 


Friction or Final Value 

Initial Value 

Inlet 

Liquid 

Node Number 

Single Phase Value 

Saturated 

Two-Phase Value 

Vapor 

Wall 
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APPENDIX B 

TRANSIENT CRYOGEN TRANSFER PROGRAM 
VERIFICATION PRESENTATION 

The presentation in this Appendix gives the results of additional 
program verification analyses which were made due to program up- 
dates since the program development reported in Letter Report 
5-9030-HT-158, "Cryogen Transfer Computer Program Development and 
Verification," dated September 3, 1974. 
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TRANSIEINIT 


CRYOGEN* 


TRANSFER 


PROGRAM 



"•'Applicable to any boiling fluid. 
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BOEING AEROSPACE COMPANY 
SPACE SYSTEMS DIVISION 
HUNTSVILLE, ALABAMA / 


CRYOGEN LOADING PROBLEM STATEMENT 


V. L. GLASGOW 
NOVEMBER 25, 1974 



$ GIVEN 

• TYPE OF FLUID (LOX, LH 2 . ETC.) 

• SYSTEM DESCRIPTION 

- PUMP PERFORMANCE OR WINDMILLING PRESSURE DROP CHARACTERISTICS 

- LINE GEOMETRY AND COMPONENTS 

- STEADY STATE HEAT LEAK 

• BOUNDARY CONDITIONS VS. TIME 

- PUMP RPM OR FLOW RATE 

- INLET CONDITIONS (PUMP INLET OR STORAGE TANK) 

- FLOW CONTROL VALVE POSITIONS 

- SYSTEM EXIT PRESSURE 


• REQUIRED 

• LOCAL TRANSIENT PROPERTIES AND FLOW CONDITIONS FOR ONE OR TWO-PHASE 
FLOW AT ANY POINT IN THE SYSTEM 


• SOLUTION 

• INTEGRATED ANALYSIS OF PUMP AND SYSTEM PERFORMANCE; FINITE ELEMENT 
METHOD INCORPORATING UNSTEADY ENERGY, MASS AND MOMENTUM WITH TWO-PHASE 
PROPERTY, FRICTION LOSS, AND HEAT TRANSFER RELATIONS 


1-06161-0810 


/ 


BOEING AEROSPACE COMPANY 
SPACE SYSTEMS DIVISION 
HUNTSVILLE, ALABAMA 


COMPUTER PROGRAM DEVELOPMENT AND VERIFICATION 


V. L. GLASGOW 
NOVEMBER 25, iy74 


\ 




TRANSIENT CRYOGEN TRANSFER 
LOADING ANALYSES 


PROGRAM (TCTP) DEVELOPED FOR SHUTTLE CRYOGEN 


DO 

I 

4 ^ 


• TCTP VERIFIED BY S-IB STAGE LOX LOADING SYSTEM ANALYSIS 

• TANK HEAD FLOW 

• PUMP RAMP FLOW 

• VERIFICATION ANALYSIS SAMPLE RESULTS 

• TANK HEAD FLOW TRANSIENTS 

• PUMP BUILD-UP FLOW AND BACK PRESSURE TRANSIENTS 

• THERMODYNAMIC STATE AS A FUNCTION OF TIME AND LOCATION 

• FLOW RATE AS A FUNCTION OF TIME AND LOCATION 

• PRESSURES AS A FUNCTION OF TIME 
AND LOCAiION 
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/ BOEING AEROSPACE COMPANY / 
^ SPACE SYSTEMS DIVISION / 
HUNTSVILLE. ALABAMA / 


S-IB LOADING SYSTEM SCHEMATIC 


S-I8 STAGE - 
LOX TANK L3 


S-IB STAGE 

ENTRANCE 
ELEVATION 71m ' 
(232.9 FT) ^ 


LOX MAST 
MEASUREMENTS 
ELEVATION 67m 
(219.7 FT) 


VALVE COMPLEX 
ELEVATION 63.5m 
(208.3 FT) 


LOX " 
STORAGE 
. TANK 


TOP OF LAUNCH PAD 
ELEVATION lO.lm 
(33.1 FT) 


ELEVATION 23.3m 
(76.6 FT)r 


0 m 
ELEV 

.(REFERENCE) 


M/L DISCONNECT 


•-310 m-H 
(1017 FT) 


-BASE OF LAUNCH PAD 
ELEVATION -0.2m 
(-0.6 FT) 


SIMULATED SYSTEM 
BOUNDARY 
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• FLOW INITIALIZED BY OPENING VALVE A28750 

• SYSTEM INITIAL TEMPERATURE DOWNSTREAM OF A28750: AMBIENT 

• SYSTEM INITIAL PRESSURE DOWNSTREAM OF A28750: 1 PSIG 

• QUALITY UPSTREAM OF A28750: 0.0 

• QUALITY DOWNSTREAM OF A28750: 1.0 

t INITIAL PUMP FLOW RATES: 0.0 

• INITIAL SYSTEM EXIT FLOW RATE: 0.0 

• AS-208 LCD DATA INPUT: PUMP INLET PRESSURE VS. TIME 

PUMP DISCHARGE TEMPERATURE VS. TIME 
SYSTEM EXIT PRESSURE VS. TIME 
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TIME FROM FLOW INITIATION 
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FROM Fin^, IMTlATIOr; 


■ F.' 
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iL-a 


f BOEING AEROSPACE COMPANY 
SPACE SYSTEMS OlVISION 
HUNTSVILLE. ALABAMA / 


S-IB PUMP RAMP TRANSIENT 


X V. L. GLASGOW 
NOVEMBER 26, 1^74 



• FIRST PUMP RPM RAMP INITIATED 636 SECONDS AFTER OPENING VALVE A28750 

# SYblEM PRESSURES, TEMPERATURES, FLOW RATES, AND QUALITIES CALCULATED 
FROM TANK HEAD FLOW 

B AS-208 LCD DATA INPUT: 

PUMP RPM VS. TIME 

PUMP INLET PRESSURE VS. TIME 

PUMP DISCHARGE TEMPERATURE VS. TIME 


SYSTEM EXIT PRESSURE VS. TIME 
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/ BOEING AEROSPACE COMPANY 
^ SPACE SYSTEMS OlVISION > 
HUNTSVILLE, ALABAMA / 


TRANSIENT SYSTEM PRESSURES 
(PUMP RAMP) 


V. L. GLASGOW 
NOVEMBER 25, 1974 


» — -L- 


^ AS-208, LCD DATA 
— ANALYTICAL DATA 


200 ■ 


I - 4- 


PUMP DISCHARGE 


M/L DISCONNECT 



^ _ 

VALVE 

COMPLEX OUTLET 

1 50 200 

250 

300 

TIME FROM INITIAL RPM 

SEC 
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VAPOR VOLUME FRACTION 


f BOEING AEROSPACE COMPANY 
SPACE SYHEMS OlViSION 
HUNTSVILLE. ALABAMA 


TRANSIENT VAPOR VOLUME FRACTIONS 
(PUMP RAMP) 



V. L. GLASGOW 
NOVEMBER 25, 197J 




ANALYTICAL VAPOR | 

VOLUME FRACTION 1 

LIQUID INDICATED BY I 
AS-208, LCD MEAS. ' 


300 


T 


ME FROM INITIAL RPw SEC 
















TEMPERATURE 


f BOEING AEROSPACE COMPANY 
SPACE SYSTEMS OlViSION 
HUNTSVILLE. ALABAMA 



SYSTEM TEMPERATURE GRADIENTS 
(PUMP RAMP) 



V. L. GLASGOW 
NOVEMBER 25, 1<?74 





D180-19190-1 


B-19 



• CURRENT COMPUTER PROGRAM CAPABILITIES 

• PREDICTS LOCAL TRANSIENT ONE OR TWO-PHASE PROPERTIES 

• ANY FLUID FOR WHICH THERMODYNAMIC TABLES ARE PROVIDED 

• VARIOUS LINE GEOMETRY, COMPONENTS AND HEAT LEAK 

• TRANSIENT TANK HEAD FLOW (WINDMILLING PUMP) 

• PUMP SPEED BUILD-UP TRANSIENTS 
t TWO-PHASE FLOW PRESSURE DROP 

• TWO-PHASE FLOW WALL HEAT TRANSFER 

• VARIABLE SYSTEM ENTRANCE PROPERTIES, PUMP PPM, AND VALVE POSITIONS 

• OPTION FOR STORAGE TANK OR PUMP INLET SYSTEM BOUNDARY 
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